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The Hohenberg-Kohn theorem of density functional theory (DFT) for the case of electrons 
interacting with an external magnetic field (that couples to spin only) is examined in more detail 
than previously. An unexpected generalization is obtained: in certain cases (which include half 
metallic ferromagnets and magnetic insulators) the ground state, and hence the spin density matrix, 
is invariant for some non-zero range of a shift in uniform magnetic field. In such cases the ground 
state energy is not a functional of the spin density matrix alone. The energy gap in an insulator or 
a half metal is shown to be a ground state property of the N-electron system in magnetic DFT. 



The half metallic state of a ferromagnet has been re- 
ceiving greatly increasing attention since its prediction 
from band theory [Q to be the ground state of impor- 
tant magnetic materials such as CrC>2, [|| NiMnSb, ||] 
and Sr2FeMoC>6 Q] and several other intermetallics and 
oxides, and its unusual physical properties. j^] Such sys- 
tems have become very attractive for magnetoelectronic 
applications, where control of the spin degree of freedom 
is already leading to new devices. j|] Materials thought to 
be half metals have been connected with the phenomena 
of colossal magnetoresistance (CMR), large tunneling 
MR, H and large, low-field intergrain MR, |J and they 
would be optimal for applications of spin valve systems 
H for non- volatile magnetic memory and for high density 
magnetic storage. 

The half metallic state is, in a one-electron picture, a 
collinear magnetic state in which one spin direction is 
metallic while the other is gapped ('insulating'). This 
state is half metallic in another sense: the absence of low 
energy spin flips leads to a vanishing magnetic suscepti- 
bility like an insulator, but its charge response (conduc- 
tivity) is that of a metal. These properties combine to 
give a one-electron description of a spin-charge separated 
state. In fact, almost all understanding of half metals so 
far is based on the one-electron picture, which opens up 
questions such as (1) what is a half metal in many body 
context, and (2) are there other unusual possibilities in 
magnetic systems? One general characterization might 
be in terms of conductivity (charge response) and sus- 
ceptibility (spin response) alluded to above: in an insu- 
lator both vanish, in a conventional (even ferromagnetic) 
metal both are non-zero, and in a half metal the con- 
ductivity is non-zero while the susceptibility vanishes. A 
clear many body formulation is however lacking. 

Since density functional theory (DFT) is a rigorous 
many body theory for (chosen) ground state properties, 
we revisit the foundations of DFT with magnetic prop- 
erties in mind. The first Hohenberg-Kohn (HK) theorem 
, which is the basis for the DFT of spin-independent 
particles, demonstrates the existence of a unique map 

n(r) i ^ v(r) mod (constant), (1) 



where v is the external potential and n is the ground state 
particle density. According to the second HK theorem, 
the ground state energy and density are obtained as the 
solution to a variational principle: 

E[v - fi\ = min n {F[n] + J n(v - fi)d 3 r}, (2) 

with [i the chemical potential. Although the variational 
principle has been put on an independent, more general 
basis, jll| the uniqueness of the map ([[]) remains an im- 
portant issue regarding the existence and uniqueness of 
the functional derivative SF/8n — —(v — fx). DFT, as ex- 
tended by Kohn and Sham [ fl2f and many others, forms 
the basis of our understanding of the electronic behavior 
of real materials. The theory has been extended to elec- 
trons with spin jl3|,[l4) and also applied heavily, however 
the HK theorem for interacting particles with spin has re- 
peatedly been stated to be analagous to the HK theorem, 
although this was already questioned in [ p"3| . Zero suscep- 
tibility, however, would imply that the ground state spin 
density does not change when an external magnetic field 
is changed. In this paper we construct a more revealing 
generalization of the HK theorem, obtain explicitly the 
conditions that allow half metallicity, and demonstrate 
some unexpected consequences. 

We consider the system in an external magnetic field 
B(r) in the (commonly considered) non-relativistic limit, 
in which the field acts only on the electron spin and the 
dipolar interaction between spins is neglected. The po- 
tentials can be combined into 2x2 spin matrix 

u ss /(r) = v(r)6 ss > - // s B(r) ■ a ss >, (3) 

The external field B may vary in magnitude and direc- 
tion. 

Realizing that two different scalar potentials cannot 
lead to the same ground state 'F, the original derivation 
of HK jl0| concluded that if 'F i— ► v — fi is unique then 
n i— » v — n is unique. Following the original derivation 
of HK pi], we begin by supposing that there are two 
different potentials u, u' which lead to the same ground 
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state Vt. We show that $ h (v — fi,B) is not a unique 
mapping in general. 

The many-body Hamiltonian of the system is 



H = T + W + U, 



(4) 



where T is the kinetic energy operator, W is the 
Coulomb interaction energy, and U is the interaction 
with the external potential. The fermionic many-particle 
Schrodinger equation is (in atomic units) 



N 



y^w(r^rj) *(riai, ..,rjvajv) 



(5) 



JV 



+ ^ ^ u ai ,Pi (r»)^(riQ!i, ..,ri/3i, ..,rjvajv) 

i ft 

= J5*(r!ai, ...,r w aiv), 

where r^, a, are the space and spin coordinates of the i-th 
electron; w(r,r ') = e 2 /|r — r'| is the Coulomb interac- 
tion. 

Assume there are two external potentials u, v! with 
energies E, E' that have the same ground state wave 
function ^(riai, rpfajsr). Subtracting the two many- 
particle Schrodinger equations leads to 

N 

'y]y^ j Au ai ^ i (T i )^(T 1 ai, ..,r f /3j, ..,r N a N ) = (6) 

i=i pi 

AE *(riai, ...,r N a N ), 

where Au = u — v! , AE = E — E' . Now, we perform a 
unitary spin rotation Q ss i (r) at each point of space that 
diagonalizes the difference in potentials (i.e. rotates B 
to lie along the z direction: 

{Q{v)[Au(v)]Q\v)} ssl = Au s (v)5 ssl . (7) 

The wavefunction is transformed according to 

N 

n^ia-fa)^ 1 "!"!' ••■' r A fQ! Jv) = *( r l«l> ■•■,r A ra A r), 



where Yi^ Qaia'.( r i) is the operator that rotates each of 
the (cti). Collecting these results gives 

N 

Ait a< (rt)#(r 1 a 1 , ...,r N a N ) = AE *(riai, r N a N ). 

i=l 

"J is some {r^j-dependent multi-component function of 
the 2 N possible spin configurations, at least one of which 
must be non-zero. Choose a non-zero component ty c and 
denote by the number of a, =t values in this compo- 
nent. Since $ is antisymmetric (as was VP) with respect to 
permutations of (r^a^) with (r^-ay), we may renumber the 



particle indices in such a way that ax = a 2 = ••• = c^jv-p 
a7v T +i = «at t +2 = ■•• = ajv- This ordering lets us write 



£ 

»=i 



N 



Au 1 (r i )+ A H(*i) >*c(ri T,-,rjv|) (8) 

i=JV T +l J 

= A£* c (n T,...,rjv i). 



This equation must hold for all values of (ri, rjy). (We 
suppose u, u' are analytic in r except possibly at isolated 
points, so that ^ c is non-zero almost everywhere.) By 
varying only r 1 , and then separately varying only rjy, we 
obtain 
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C T , Am 



C 



1 ! 



(9) 



where C-f , Cy are constants. The special cases N-[ = or 
N-[ = N do not lead to new consequences. For further 
analysis, we consider separate cases. 
Case A: impure spin states. Suppose that there are at 
least two components of ^> with different values of iV-j- 
and hence N l = N - Np Then 



iV T C* T + (JV - N^Ct = AE. 



(10) 



Since this holds for two different values of Nf, it fol- 
lows that Cf = C \ = C, which leads to Aw-f = ASj so 
Au a ^ = C5 a< 0. By the ground state energy minimum 
principle, this recovers the usual Hohenberg-Kohn result 



v(r) - v'(r) = C, 
B(r) -B'(r) = 0, 



(11) 



implying a non-zero susceptibility. 

Case B: pure spin states. Suppose now that all non-zero 
components of ^ have the same value of and N±. 
These may be considered as "pure spin" states, eigen- 
functions of the operator S z = J^i a ai8i/^ w ith eigenval- 
ues S z = Nf - (N/2) = (JV T - Ni)/2. 'Then C* T and Cj. 
need not be equal and we can write 



Ait 



C T 

o c l 



CI- fi B Bcr z , 



(12) 



where C = (C* T + C\)/2 and -[l b B = (C T - C x )/2. 

Backtransforming according to the inverse of Eq. (Q) 
gives 



Au a p(v) = CSa/3 + hbB[Q^ (r)a z Q(r)] a p. 



(13) 



The last term on the right is position-dependent, non- 
diagonal, and non-vanishing in general. In this case the 
conditions for identical ground state wavefunctions are 



v(r)-v'(r) = C 
B(r) - B'(r) = B e(r) 



(14) 



where e is the unit vector ^Tr{aQ^(r)a z Q(r)}. The re- 
sult (11) is modified accordingly. This result is a highly 
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non-trivial generalization of the HK theorem: two mag- 
netic fields whose difference is constant in magnitude, but 
possibly is non-unidirectional, may give rise to the same 
ground state. 

Now we investigate the conditions on Q for which ^> is 
an eigenstate of S z , i.e. $ describes a collinear spin ar- 
rangement. Considering the Hamiltonian Eq. (^), there 
must be an operator 



U 



N 

E« 

i=l 



Ta'.CH ( r i) a a'.0'.Q/3'M T i) 



(15) 



that commutes with T + U + W. We now specialize to 
the particular case where one of the external fields, B', is 
zero. Since the interaction W is spin independent, U will 
commute with H' if and only if it commutes with T. One 
can show that this necessitates that Q be r- independent, 
so that VP itself is an eigenstate of S z , and hence is a 
collinear spin state. The second condition in Eq. Jl4| ) 
reduces to B — B' = Bz: a turning on of a uniform mag- 
netic field leaves the ground state invariant. Restated: 
in the subspace of collinear magnetizations, the ground 
state determines the magnetic field only up to some codi- 
rectional uniform field. A direct corollary is that there is 
no longer any ground state energy functional E[n] of the 
density n ss i (r) alone (see Eq. ( jig ) below) . 

In the collinear situation, inserting Au = [1bBo~ z into 
Eq. (6) yields 



AE(B) = -(Ni-N l )n B B ) 



(16) 



which gives the well known dependence of energy vs. field 
for a system of fixed spin. Consider as a simple example 
a Be atom in a uniform magnetic field, with its ground 
state characterized as ls 2 2s 2 (N = 4, Nf = 2). The 
lowest excited ^-eigenstate is ls 2 2s2p with iVj = 3. Its 
excitation energy is that of a 2s — > 2p promotion. There 
is another excited state ls2s 2 2p with the same N^ t but 
the much higher excitation energy of a Is — > 2p core exci- 
tation. The energetically lowest — 4 state is Is2s2p 2 
whose excitation energy is roughly the sum of the previ- 
ous two. The situation is sketched in Fig. 1(a), where 
the lines with positive slopes correspond to states with 
all spins reversed. 

Since states with JV-j- = N/2 i n are degenerate for 
B = 0, Fig. 1(a) may be supplemented symmetrically to 
the vertical axis. Hence, for \B\ < Bq the groundstate 
is ls 2 2s 2 with energy Eq, for Bo < B < B\ the ground 
state is ls 2 2s2p with energy E\ — 2/ibB, and for B > B\ 
the ground state is \s2s2p 2 with energy £3 — 4/is_B. The 
ground state does not change with field except at certain 
isolated values. 

In an extended system, say a non-magnetic insulator 
with gap A g , there is a continuum above A g (one excited 
electron with reversed spin), another continuum above 
2A g (two excited electrons) and so on, as illustrated in 



Fig. 1(b). In an extended system one would prefer to 
consider the intensive quantity 



AE(B) fty-Ni 



N 



N 



(17) 



instead of AE itself. Then, one finds that for hbB < A g 
the groundstate is independent of B, beyond which the 
state changes and AE /N veers off. Thus while the gap 
A 9 is not a ground state property of the N particle sys- 
tem in paramagnetic DFT (it involves the Nil particle 
ground states), it is a ground state property in the pres- 
ence of a uniform field. 

For a stoichiometric half metal with moment per cell 
Hb-M. (Ai an integer) the picture is related, except there 
is an overall bias - a slope of -\xbM. in the energy per 
cell - and the positive and negative B directions are not 
symmetric. The situation that is sketched in Fig. 1(c) 
has a gap A„ + A c for J, spin states, with no gap for 
f spin. The chemical potential fi corresponds to the en- 
ergy to remove an f spin, and the quantities A t , = fisB v , 
A c = hb\B c \ represent the energy, or field, required to 
flip a spin from | to |, or vice versa. Note again that the 
interval of B for which the state does not change, which 
is the gap in the j spectrum, is a ground state property 
of the N particle system in an external magnetic field. 

It is useful to consider the form of constrained DFT in 
which Nij and N± are specified, which leads to two asso- 
ciated chemical potentials /xj, /ij. Then as N s is changed 
to N s i 1, fjb s may vary only to order 1/N S (metallic be- 
havior) or it may jump discontinuously across a gap, just 
as is the case for insulators. jl5| The half metal is de- 
fined as that situation in which one and only one of fi s 
(we have choosen |) is discontinuous upon addition of 
one electron. For an insulator, there is a discontinuity in 
/i for both spins. 

We now consider the KS eigenvalue spectrum. As long 
as the external field shifts the bands sufficiently little 
not to disturb the half metallicity (B c < B < B v ), 
the ground state, and hence the charge density in each 
spin channel, remain unchanged. Using the same argu- 
ments as were applied to establish the discontinuity in 
Vxc(N(/j>)) for an insulator as /i crosses the gap (the ki- 
netic energy is discontinuous across the gap) [ jUj , one 
finds that there is a discontinuity in v X c. i (Ni, N±) if the 
filling with N± moves (i± across the gap. |l6| 

The Kohn-Sham gap s g i is smaller than the true 
(quasiparticle) gap A g = A c + A v . When the magnetic 
field is large enough that fi reaches the KS band min- 
imum e c (Ni), the occupation of that channel becomes 
Nj, + e (with e — > 0). This is the point of the disconti- 
nuity, where the KS conduction eigenvalue (in fact, the 
entire j spectrum) jumps upward. By comparison with 
Dyson's equation, and the fact that the system's ground 
state spin densities must be the same whether obtained 
from DFT or the quasiparticle Greens function, this jump 
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must be such as to make e c i(N^, iVj_ + e) = A c , the quasi- 
particle conduction band edge, for e — >0. 

It is apparent then that the KS gap in the insulating 
channel is not equal to the true gap in that channel, and 
that e c (A^-|-, iVjJ — //is not the true spin flip energy (which 
is A c - /i). By our definition, |u| as the reverse field is 
applied and fi is driven toward the valence band maxi- 
mum Svi, there is no discontinuity, and the other spin 
flip energy - a true excitation energy - is given correctly 
by DFT. Needless to say, an approximation such as the 
local density approximation that interpolates across the 
discontinuity, will fail to predict both A c and A„. 

We now summarize. We have presented new, rigorous 
results for the Hohenberg-Kohn mapping in a magnetic 
field. We obtain conditions that characterize half metals: 
(1) two collinear systems in different uniform magnetic 
fields may have the same half metallic (or magnetic in- 
sulating) ground state; (2) exactly one of the chemical 
potentials /i s is discontinuous upon particle addition to 
a half metal. We have pointed out other consequences, 
primary among them being that the ground state en- 
ergy of a system is no longer a unique functional of the 
density n S8 t when magnetic fields are allowed (although 
the ground state itself is), and that the gap in a half 
metal is a ground state property of the N particle sys- 
tem. These results are only exact in the non-relativistic 
(c — > oo) limit. For c finite, half metallicity is an approxi- 
mate notion due to orbital currents and orbital moments 
and spin-orbit coupling that mixes them, and the general 
theory jl7| probably restores the conventional theorems 
of DFT. Still, the notion of half metallicity will be an 
important model limit. 
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FIG. 1. Energy change with field B for: (a) Be atom. The levels are discrete and the variation is —S z ^bB. (b) a nonmagnetic 
insulator. The ground state energy is constant over a range of fields — A g < hbB < A g , beyond which magnetization is induced, 
(c) a half metal. The ground state energy varies as — /Us(iVf — N±)B in the range B c < B < B v (A v < hbB < A c ). In (b) and 
(c), the small dots indicate a continuum extending upwards. See text for further explanation. 
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